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We study non-local two-qubit operations from a geometric perspective. By applying a Cartan 
decomposition to su(4), we find that the geometric structure of non-local gates is a 3-Torus. We 
derive the invariants for local transformations, and connect these local invariants to the coordinates 
of the 3- Torus. Since different points on the 3-Torus may correspond to the same local equivalence 
class, we use the Weyl group theory to reduce the symmetry. We show that the local equivalence 
classes of two-qubit gates are in one-to-one correspondence with the points in a tetrahedron except 
on the base. We then study the properties of perfect entanglers, that is, the two-qubit operations 
that can generate maximally entangled states from some initially separable states. We provide 
criteria to determine whether a given two-qubit gate is a perfect entangler and establish a geometric 
description of perfect entanglers by making use of the tetrahedral representation of non-local gates. 
We find that exactly half the non-local gates are perfect entanglers. We also investigate the non- 
local operations generated by a given Hamiltonian. We first study the gates that can be directly 
generated by a Hamiltonian. Then we explicitly construct a quantum circuit that contains at most 
three non-local gates generated by a two-body interaction Hamiltonian, together with at most four 
local gates generated by single qubit terms. We prove that such a quantum circuit can simulate 
any arbitrary two-qubit gate exactly, and hence it provides an efficient implementation of universal 
quantum computation and simulation. 



I. INTRODUCTION 



Considerable effort has been made on the characterization of non-local properties of quantum states and operations. 
Grassl et al. [jj] have computed locally invariant polynomial functions of density matrix elements. Makhlin has 
recently analyzed non-local properties of two-qubit gates and presented local invariants for an operation M G U(4). 
Makhlin also studied some basic properties of perfect entanglers, which are defined as the unitary operations that 
can generate maximal entangled states from some initially separable states. Also shown were entangling properties of 
gates generated by several different Hamiltonian operators. All these results are crucial for physical implementations 
of quantum computation schemes. 

Determining the entangling capabilities of operations generated by a given physical system is another intriguing and 
complementary issue. Zanardi || has explored the entangling power of quantum evolutions. The most extensive 
recent effort to characterize entangling operations is due to Cirac and coworkers ||, [g, [?], |[ [| |o[ O, [l2| . Kraus and 
Cirac Q focused on finding the best separable two-qubit input states such that some given unitary transformation 
can create maximal entanglement. Vidal et al. ^ developed the interaction cost for a non-local operation as the 
optimal time to generate it from a given Hamiltonian. The same group, Hammerer et al. |12| then extended these 
considerations to characterize non-local gates. These works are closely related to time optimal control as addressed 
recently by Khaneja et al. [T^ |, who studied systems described by a Hamiltonian that contains both a non-local internal 
or drift term, and a local control term. All these studies assume that any single qubit operation can be achieved 
almost instantaneously. This is a good approximation for the situation when the control terms in the Hamiltonian 
can be made large compared to the internal couplings. 

Universality and controllability are issues of crucial importance in physical implementations of quantum information 
processing 
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series of important results have been obtained since questions of universality were first addressed 
by Dcutsch in his seminal papers on quantum computing j^J. Deutsch ]l7[ | proved that any unitary operation 
can be constructed from generalized Toffoli gate operating on three qubits. DiVincenzo |l8| proved universality for 
two-qubit gates by reconstructing three-qubit operations using these gates and a local NOT gate. Similarly, Barenco 
fl9| , |20|| and Sleator and Weinfurter identified the controlled unitary operation as a universal two-qubit gate. 
Barenco |^2j showed the universality of the CNOT gate supplemented with any single qubit unitaries, and pointed 
out advantages of CNOT in the context of quantum information processing. Lloyd p3[ showed that almost any 
quantum gate for two or more qubits is universal. Deutsch et al. p4j proved that almost any two-qubit gate is 
universal by showing that the set of non-universal operations in C/(4j is of lower dimension than the U(4) group. 
Universal properties of quantum gates acting on an n > 2 dimensional Hilbert space have been studied by Brylinski 
p5| . Dodd et al. p(| have pointed out that universal quantum computation can be achieved by any entangling 
gate supplemented with local operations. Bremner et al. |^7|] recently demonstrated this by extending the results 



2 



of Brylinski, giving a constructive proof that any two-qubit entangling gate can generate CNOT if arbitrary single 
qubit operations are also available. Universal sets of quantum gates for n-qubit systems have been explored by Vlasov 
p8| . in connection with Clifford algebras. 

General results on efficient simulation of any unitary operation in SU(2 n ) by a discrete set of gates are embodied 
in the Solovay-Kitaev theorem ^o) and in recent work due to Harrow et al. f3lfl . The Solovay-Kitaev theorem 
implies the equivalence of different designs of universal quantum computers based on suitable discrete sets of single- 
qubit and two-qubit operations in a quantum circuit. An example is the standard universal set of gates including 
controlled-NOT and three discrete single-qubit gates, namely Hadamard, Phase and 7r/2 gates p5| . Other universal 
sets have also been proposed. According to the Solovay-Kitaev theorem, every such design can represent a circuit 
that is formulated using the standard set of gates. Consequently, all quantum computation constructions — including 
algorithms, error-correction, and fault-tolerance — can be efficiently simulated by physical systems that can provide a 
suitable set of operations, and do not necessarily need to be implemented by the standard gates. This moves the focus 
from study of gates to study of the Hamiltonians whose time evolution gives rise to the gates. In this context, Burkard 
et al. |]32f studied the quantum computation potential of the isotropic exchange Hamiltonian. This interaction can 
generate vSWAP gate directly. However, CNOT cannot be ob tained directly from the exchange interaction. Burkard 
et al. showed that it can be generated via a circuit of two VSWAP gates and a single-qubit phase rotation. More 
recently, Whaley and co-workers have shown that the two-particle exchange interaction is universal when physical 
qubits arc encoded into logical qubits, allowing a universal gate set to be constructed from this interaction alone 
p3| , 34, |35|, [36| |37|, |38|]. This has given rise to the notion of "encoded universality", in which a convenient physical 



interaction is made universal by encoding into a subspace [p5[ |3q] . Isotropic, anisotropic, and generalized forms of the 
exchange interaction have recently been shown to possess considerable power for efficient construction of universal 
gate sets, allowing explicit universal gate constructions that require only a small number of physical operations 

In this paper, we analyze non-local two-qubit operations from a geometric perspective and show that considerable 
insight can be achieved with this approach. We are concerned with three main questions here. First, achieving a 
geometric representation of two-qubit gates. Second, characterizing or quantifying all operations that can generate 
maximal entanglement. Third, exact simulation of any arbitrary two-qubit gates from a given two-body physical 
interaction together with single qubit gates. The fundamental mathematical techniques we employ are Cartan decom- 
position and Weyl group in the Lie group representation theory. The application of these theories to the Lie algebra 
su(4) provides us with a natural and intuitive geometric approach to investigate the properties of non-local two-qubit 
operations. This geometric approach reveals the nature of the problems intrinsically and allows a general formulation 
of solutions to the three issues of interest here. 

Because it is the non-local properties that generate entanglement in quantum systems, we first study the invariants 
and geometric representation of non-local two-qubit operations. A pair of two-qubit operations are called locally 
equivalent if they differ only by local operations. We apply the Cartan decomposition theorem to su(4), the Lie 
algebra of the special unitary group SU(4). We find that the geometric structure of non-local gates is a 3- Torus. 
On the other hand, the Cartan decomposition of su(4) derived from the complexification of sl(4) yields an easy way 
to derive the invariants for local transformations. These invariants can be used to determine whether two gates are 
locally equivalent. Moreover, we establish the relation of these local invariants to the coordinates of the 3- Torus. 
This provides with a way to compute the corresponding points on the 3-Torus for a given gate. It turns out that a 
single non-local gate may correspond to finitely many different points on the 3- Torus. If we represent these points 
in a cube with side length it, there is obvious symmetry between these points. We then use the Weyl group theory 
to reduce this symmetry. We know that in this case the Weyl group is generated by a set of reflections in R 3 . It is 
these reflections that create the kaleidoscopic symmetry of points that correspond to the same non-local gate in the 
cube. We can explicitly compute these reflections, and thereby show that the local equivalence classes of two-qubit 
gates are in one-to-one correspondence with the points in a tetrahedron except on the base. This provides a complete 
geometric representation of non-local two-qubit operations. 

The second objective of this paper is to explore the properties of perfect entanglers, that is, the quantum gates that 
can generate maximally entangled states from some initially separable states. We start with criteria to determine 
whether a given two-qubit gate is a perfect entangler. A condition for such a gate has been stated in We 
provide here a proof of this condition and show that the condition can be employed within our geometric analysis to 
determine which fraction of all non-local two-qubit gates are perfect entanglers. We show that the entangling property 
of a quantum gate is only determined by its geometric representation on the 3-Torus. Using the result that every point 
on the tetrahedron corresponds to a local equivalence class, we then show that the set of all perfect entanglers is a 
polyhedron with seven faces and possessing a volume equal to exactly half that of the tetrahedron. This implies that 
amongst all the non-local two-qubit operations, exactly half of them are capable of generating maximal entanglement. 

Finally, we explore universality and controllability aspects of non-local properties of given physical interactions 
and the potential of such specified Hamiltonians to generate perfect entanglers. Our motivation is related to that of 



FIG. 1: Partition of all the gates in SU(4) 



encoded universality, namely, determining the potential for universal quantum computation and simulation of a given 
physical Hamiltonian. However, whereas encoded universality sought to construct encodings to achieve universality 
of quantum logic, here we focus on the simulation of any arbitrary two-qubit gate. Achievement of this, together 
with our second result above, allows generation of maximal entanglement as well as providing universality. To realize 
this, we consider here the conventional scenario of a Hamiltonian acting on a physical set of qubits such that any 
arbitrary single qubit operation and certain specific two-qubit operations may be turned on for selected time durations 
in series. Generally speaking, two-qubit interactions include both local and non-local terms. The non-local terms can 
give rise to not only well-known entangling gates such as CNOT, but also to many other classes of gates that may 
or may not lie in the perfect entangling sector. We therefore seek a systematic way to construct quantum circuits 
from a given physical Hamiltonian that can simulate any arbitrary two-qubit gate exactly, without restriction by 
limitations of speed due to time-scales of control interactions (e.g. as in NMR) or to fundamental considerations [|40| . 
As in the study of encoded universality we start with the gates that can be directly generated by a given Hamiltonian. 
Generally, these gates form a one dimensional subset on the 3- Torus. To construct an exact simulation of any arbitrary 
two-qubit gate, we make use of the quantum circuit model. We explicitly construct a quantum circuit that contains 
three non-local gates generated by a given two-body interaction Hamiltonian for corresponding finite time durations, 
together with at most four local gates. We prove that such a quantum circuit can simulate any arbitrary two-qubit 
operation exactly and is therefore universal. In particular, it can therefore efficiently provide maximal entanglement 
from any arbitrary Hamiltonian of this form. Such efficient construction from any given Hamiltonian is extremely 
useful for design and experimental implementation of quantum information processing schemes. 



II. PRELIMINARIES 



In this section, we briefly review some basic facts about Cartan decomposition and the Weyl group within Lie group 
representation theory ^l], , and then apply these results to su(4) , the Lie algebra of the special unitary group 

SU (4). Applications of Cartan decomposition to quantum system control can also be found in |L3 |. 

We concentrate on SU(A) when studying two-qubit gates. It is well-known that an arbitrary two-qubit gate 
Ua € U(4) can be decomposed as the product of a gate U\ G SU(4) and a global phase shift e la , where a S M. 
Because the global phase has no significance in quantum mechanics, we can thereby reduce the study of the group 
C/(4) of two-qubit quantum evolution operators to S£/(4). Extensions of results from the group ST/ (4) back to Z7(4) 
are made when appropriate. 

We heuristically introduce a partition of the set of two-qubit operations represented by the group St/(4). This set 
splits into two subsets, one of local gates SU(2) ® SU(2) and the other of non-local gates SU(4)\SU(2) ® SU(2). 
The latter splits further into a set of perfect entanglers, i.e., those that can generate maximally entangled states, an 
example of which is CNOT, and the complementary set of those non-local gates that are not perfect entanglers. This 



schematic partition is illustrated in Fig. [H. A rigorous definition of perfect entanglers is presented in Sec. 



IV. 
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A. Cartan decomposition and Weyl group 



Our first goal is to establish fundamentals for a geometric picture of non-local unitary operations with emphasis 
on their generators, which are represented by Hamiltonian operators in physical context. We start with a summary 
of some basic definitions |^|, |43). Consider a Lie group G and its corresponding Lie algebra g. The adjoint 
representation Ad g is a map from the Lie algebra g to g which is the differential of the conjugation map a g from 
the Lie group G to G given by a g (h) = ghg^ 1 . For matrix Lie algebras, Ad g (Y) = gYg^ 1 , where g, Y are both 
represented as matrices of compatible dimensions. The differential of the adjoint representation is denoted by ad, and 
a.dx is a map from the Lie algebra g to g given by the Lie bracket with X, that is, adxiY) = [X, Y]. 

We now define an inner product on g by the Killing form B(X, Y) = tr(adx ady). Let {X\, . . . , X n } be a basis for 
such that 



The numbers Cj fc <E 



[Xj,Xk] — ^ C l jk Xi 



(1) 



are the structure constants of the Lie algebra g with respect to the basis, where j, k run from 1 to n. Since 

n n 

adx 3 [X x , ■ • • , X n ] = Cj t Xi, . . . , ^2 Cjn-X-i] 



1=1 



• X n ] 



=i 
i 



r 1 ' 

jn 



(2) 



the matrix representation of adj^. with respect to the basis is 



r 1 ' 



f~vn 



(3) 



Thus, the trace of adxj &&x k , which is B(X 3 ■, Xk), is J2ab=i Cj b Cl a , which is also the jfc-th entry of the matrix of 
the quadratic form B{ , ). The Lie algebra g is semisimple if and only if the Killing form is nondegenerate, i.e., the 
determinant of its matrix is nonzero. 

Let A' be a compact subgroup of G, and t the Lie algebra of K . Assume that g admits a direct sum decomposition 
g = p © t, such that p = 1 with respect to the metric induced by the inner product. 

Definition 1 (Cartan decomposition of the Lie algebra g) Let g be a semisimple Lie algebra and let the de- 
composition g = p © t, p = t satisfy the commutation relations 

[e,e]ce, [p,e]cp, [p,p]ct (4) 

This decomposition is called a Cartan decomposition of g, and the pair (g, t) is called an orthogonal symmetric Lie 
algebra pair. 

A maximal Abclian subalgebra o contained in p is called a Cartan subalgebra of the pair (g, t). If a' is another Cartan 
subalgebra of (g,6), then there exists an element k £ K such that Adfe(o) = a'. Moreover, we have p = U^eK Adfe(a). 

Proposition 1 (Decomposition of the Lie group G) Given a semisimple Lie algebra g and its Cartan decompo- 
sition g = p © I, let a be a Cartan subalgebra of the pair (g, 6), then G = A"exp(o)A. 



For X £ a, let W £ g be an eigenvector of adx and a(X) the corresponding eigenvalue, i.e., 

[X, W] = a{X)W. 



(5) 



The linear function a is called a root of g with respect to o. Let A denote the set of nonzero roots, and A p denote 
the set of roots in A which do not vanish identically on a. Note that if a € A, it is also true that —a £ A. 
Let M and M' denote the centralizer and normalizer of a in K, respectively. In other words, 



M = {k£ K\ Ad fe (A) = X for each X £ a}, 
M' = {k£ K\ Ad fc (a) c a}. 



(6) 



5 



Definition 2 (Weyl group) The quotient group M'/M is called the Weyl group of the pair (G,K). It is denoted 
byW(G,K). 

One can prove that W(G, K) is a finite group. Each a <G A p defines a hyperplane a(X) = in the vector space a. 
These hyperplanes divide the space a into finitely many connected components, called the Weyl chambers. For each 
a 6 A p , let s a denote the reflection with respect to the hyperplane a(X) — in a. 

Proposition 2 (Generation of the Weyl group) The Weyl group is generated by the reflections s a , a £ A p . 
This proposition is proved in Corollary 2.13, Ch. VII in | fH|| . 

B. Application to su(4) 

Now we apply the above results to su(4), the Lie algebra of the special unitary group SU(4). The Lie algebra 
= su(4) has a direct sum decomposition g = p ffi t, where 

« = span Ual, a\, a 2 x , a 2 y , a 2 }, 

(7) 

^ r 1 2 12 12 12 12 12 12 12 12t 

p = span -{a x a x , a x a y , a x a z , a y a xl a y a yl a y a zl a z a x , a z a y , a z a z ). 

Here a x , o~ y , and o~ z are the Pauli matrices, and cr^cri = <r\ ® o~\. If we use Xj to denote the matrices in Eq. (]?]), 
where j runs from left to right in Eq. (j?j), we can derive the Lie brackets of Xj and Xk- These are summarized in 
the following table: 



[Xj, Xk] 


x x 


x 2 


x 3 


x 4 


x 5 


x 6 


x 7 


X 8 


x 9 


^10 


Xn 


-Xl2 


X\ 3 


Xxi 


■X~15 


X x 





-x 3 


x 2 




















— Xx 3 


—Xu 


— X\§ 


Xx 9 


Xn 




x 2 


X 3 





-X x 











-^13 


Xu 


*lfi 











-x 7 


-X 8 


-x 9 


x 3 


-x 2 


Xi 














— Xio 


-Xn 


—^12 


x 7 


X S 













x 4 














-x 6 


x 5 





—X 9 


X 8 





— X\2 







—-^15 


Xxi 


x 5 











x 6 





-Xi 


x 9 





-x 7 


X\ 2 





— X\Q 


Xx§ 





— Xi 3 


x e 











-x 5 


Xi 





-X 8 


x 7 





-Xn 


XlO 





—X\i 


-^13 





x 7 





— X\3 


Xxo 





-Xg 


^8 





-X 6 


*5 


-x 3 








x 2 








X s 





— X14 


Xn 


Xg 





-x 7 







-x. 





-x 3 








x 2 





x 9 





— -X15 


X\ 2 


-X s 


x 7 





-x 5 


Xi 











-*3 








x 2 


X\o 


X\3 





-x 7 





—X%2 


Xn 


x 3 











-x 6 




-Xx 








X n 


Xu 





-X 8 


Xi 2 





—Xio 





x 3 





x 6 





-Xi 





-Xx 





X\2 


X15 





-x 9 


-Xn 


^10 











x 3 


-x 5 


Xi 











-Xx 


X\z 


—Xio 


x 7 








— Xis 


Xu 


-x 2 








Xx 











-x 6 


x 5 


X\i 


-X n 


x s 





x 15 





—Xi 3 


















x 6 





-Xi 


X\§ 


—X%2 


x 9 





—Xn 


x 13 











-x 2 








X x 


-x 5 


x 4 






Now the structure constants Cj fe can be found from the above table (see Eq. ([[])) so that we can evaluate 



15 15 

B(Xj, X k ) =J2J2 C !" C L = -8S jk . (8) 

a=l 6=1 

It is easy to verify that tr(XjXk) — —fijk, and thus the Killing form of su(4) is B(X,Y) = 8tr(XY). Since 6 = 
span{Xi, . . . ,X§\ and p = spanjX?, . . . ,Xls}, from the Lie bracket computation table above, it is clear that 

[Met, [p,e]c P , [p,p]ce. (9) 

Therefore the decomposition g = 6© p is a Cartan decomposition of su(4). Note that the Abelian subalgebra 



a = span -{cr>^, cr^, cr^} 



(10) 
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is contained in p and is a maximal Abelian subalgebra, i.e., we cannot find any other Abelian subalgebra of p that 
contains a. Hence it is a Cartan subalgebra of the pair (fl, £). Further, since the set of all the local gates K is a 
connected Lie subgroup SU(2) £g> SU(2) of SU(A), and there is one-to-one correspondence between connected Lie 
subgroups of a Lie group and subalgebras of its Lie algebra 0], it is clear that t in Eq. (R) is just the Lie subalgebra 
corresponding to K. From Proposition [l|, any U G SU(4) can be decomposed as 

U = k x Ak 2 = fa exp{^(cicr^ + c 2 o-*o^ + c 3 cr]crf )}k 2 , (11) 

where fc 2 € S77(2) ® 51/(2), and Ci, c 2 , c 3 € M. 

Another more intuitive Cartan decomposition of su(4) can be obtained via the complexification of sl(4). Consider 
G = SX(4), the real special linear group, and A" = 50(4), the special orthogonal group. The Lie algebra st(4) is 
the set of 4 x 4 real matrices of trace zero, and so (4) the set of 4 x 4 real skew symmetric matrices. Then st(4) can 
be decomposed as st(4) = so(4) © p, where p is the set of 4 x 4 real symmetric matrices. This is nothing but the 
decomposition of a matrix into symmetric and skew symmetric parts, and it is indeed a Cartan decomposition of 
st(4). Consider the following subset of the complexification of sl(4): 

0M = so(4) + ip. (12) 

It can be verified that g M is exactly su(4), and thus (g M ,so(4)) is an orthogonal symmetric Lie algebra pair. The 
isomorphism carrying t in Eq. (|7|) into so (4) is just the transformation from the standard computational basis of 
states to the Bell basis in p2l . This procedure is of crucial importance in computing the invariants for two-qubit 



gates under local transformations. See Section |I| for more details. 

Now let us compute the Weyl group W(G, K). Let X = |(cict^ct^ + c 2 o-y<jy + c-io\a 2 z ) <E a. Identify a with R 3 , then 
X = [ci,c 2 , C3]. The roots of g with respect to a are eigenvalues of the matrix of ad^: 

A p = i{ci - c 2 , -C\ ~ c>x,—c\ - c 3 , a - c 3 , c 2 - c 3 , c 2 + c 3 , 



ci + c 2 , c\ + c 2 , ci + c 3 , -ci + c 3 , -c 2 + c 3 , -c 2 - c 3 }. 



(13) 



For a = i(ci — c 3 ) S A p , the plane a(X) = in a is the set {A G M 3 |m t A = 0}, where u = [1, 0, — 1] T . The reflection 
of X = [ci, c 2 , c 3 ] with respect to the plane a(X) — is 

2uu T 

s a (X)=X- in ^X=[c 3 ,C2,c 1 }. (14) 

Similarly, we can compute all the reflections s a as follows: 

s i{c 3 -c 2 )(X) = [ci,c 3 ,c 2 ], Sj( C2+C3 )(X) = [ci,-c 3 ,-c 2 ], 

S i(c 2 -ci) W = [C2,C1,C 3 ], S i(ci+C3 )(X) = [-C 2 ,-Ci,C 3 ], (15) 

Si( Cl _ C3 )(X) = [c 3 ,c 2 ,ci], s,( Cl+C3 )(A) = [-c 3 ,c 2 , -ci]. 

From Proposition ||, the Weyl group W(G, K) is generated by s a given in Eq. (|l5|). Therefore, the reflections s Q are 
equivalent to either permutations of the elements of [ci, c 2 , c 3 ], or permutations with sign flips of two elements. 



III. NON-LOCAL OPERATIONS 



We now study non-local two-qubit operations within the group theoretical framework of the previous section. The 
Cartan decomposition of su(4) provides us with a good starting point to explore the invariants under local gate 
operations. It also reveals that the geometric structure of the local equivalence classes is none other than a 3- Torus. 
Every point on this 3- Torus corresponds to a local equivalence class of two-qubit gates. Different points may also 
correspond to the same equivalence class. To reduce this symmetry, we apply the Weyl group theory. We show that 
the local equivalence classes of two-qubit gates are in one-to-one correspondence with the points in a tetrahedron, 
except on the base where there are two equivalent areas. This tetrahedral representation of non-local operations plays 
a central role in our subsequent discussion of perfect entanglers and the design of universal quantum circuits. 



A. Local invariants and local equivalence classes 



Two unitary transformations U, U\ S SU{4) are called locally equivalent if they differ only by local operations: 
U = k\U\k 2 , where ki, k 2 S SU(2) <§5 SU(2) are local gates. This clearly defines an equivalence relation on the Lie 
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group SU(4). We denote the equivalence class of a unitary transformation U as [U]. From the Cartan decomposition 



of su(4) in Section II B , any two-qubit gate U £ SU (4) can be written in the following form 



U = k 1 Ak 2 = fa exp{- (aerial + 020^0^ + c 3 ala 2 )}k 2 , (16) 

where kx, ki £ SU(2) eg) SU(2). Because the two-qubit gate U is periodic in Ck and the minimum positive period is 
7r, the geometric structure of \c\, C2, C3] is a 3-Torus, T 3 = S 1 x S 1 x S 1 . 

In M, local invariants were given for two-qubit gates. Here we will connect these invariants of Makhlin to the 
coordinates [01,02,03] on the 3- Torus. We first consider the case of the two-qubit gates in 5£/(4), and then extend 
the results to the general case of U(A). 



SU(4) Operations 



Consider the transformation from the standard basis of states 1 00) , 1 01) , 1 10) , 1 1 1) to the Bell basis |$ + ) = 
^(|00) + |11)), |$-) = ^(|01) + jlO)), |*+) = ^(|01) - |10», I*-) = ^(|00) - |11)). In this basis, the two-qubit 



Q 'T2 



(18) 



gate U in Eq. (16) can be written as 

U B = QtUQ = Q ] faAk 2 Q, (17) 

where 

(l i \ 
Oil 
i -1 
\1 -%) 

Recalling that |{ct^., crl, cr*, er 2 , a 2 , a 2 } is a basis for t, it is not hard to verify that |Q^{o"^, oi, <j\, a 2 , a 2 , a 2 }Q forms 
a basis for so (4), the Lie algebra of the special orthogonal group 50(4). Hence Ub can be written as 

U B = 1 Q^AQ0 2 , 

where 

01 = Q)faQ £ 50(4), 

02 = Q i k 2 Q e 50(4). 

Eq. ( p^| ) can also be obta ined from the Cartan decomposition of su(4) derived from the complexification of sl(4), as 
discussed in Section [IB. An Abelian subalgebra o is generated by \{a\a'l.,a\ j a 2 ,a\a 2 z \, and the transformation to 
the Bell basis takes these operators to §{0*, —of, ^l^ 2 }- Therefore, we have Ub = O1FO2, where 

F = Q^AQ 



(19) 



(20) 



diagje 2 j e l 2 



c l+ c 2+ c 3 



}• 



Let 



where O2 is defined by Eq. 
spectrum of the matrix m 



= UWb = Oi F 2 



(21) 



(22) 



1). The complete set of local invariants of a two-qubit gate U £ 5C/(4) is given by the 
and hence by the eigenvalues of F 2 : 



{■ 



t(ci — C2+C3) „i(ci+C2— C3) p — i(ci+C2+C3) i{ — C1+C2+C3) 



}• 



Since m is unitary and det m = 1, the characteristic polynomial of m is then 

1 



\sl — m\ = s 4 — tr(m) s 3 + — (tr 2 (™) — tr(™ 2 ))s 2 — tr(m) s + 1. 



(23) 



(24) 



Therefore the spectrum of m is completely determined by only the two quantities tr(m) and tr 2 (m) — tr (m 2 ). For a 
two-qubit gate U given in Eq. (lq), its local invariants can be derived from Eq. (G3h as: 



tr (m) =4 cos c\ cos C2 cos C3 + 4z sin ci sin 02 sin C3 , 
tr (m) — tr (m ) = 16 cos c\ cos 02 cos C3 — 16 sin c\ sin C2 sin C3 — 4 cos 2ci cos 2c2 cos 2C3 



(25) 
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2. Generalization to U(4) 

Now let us consider the local invariants for the general case of f/(4) gj. An arbitrary two-qubit gate [7 € C/(4) can 
be decomposed as the product of a gate E/i £ ££7(4) and a global phase shift e la , where det U — e l4a . It follows that 
miJJx) = e- l2a m(U), whe 



re 



m(U) = {Q^UQY QWQ, (26) 

and 

tr(m(U 1 ))=e-* 2a tr(m(U)), 
tr 2 (m(C/i)) - tr (m 2 (f/i)) = e- 44Q (tr 2 (m([/)) - tr(rn 2 ([/))). 

It is clear that the global phase factor just rotates the eigenvalues of m(U) along the unit circle in the complex 
plane, while keeping their relative phase invariant. Therefore, it does not affect the entangling properties and we can 
consequently divide by det(f/). The local invariants of a two-qubit gate U are thus given by 



Gi - 
G 2 = 



tr 2 (m{U)) 

16 det U ' (2g) 
tr 2 (m(C/)) -tr(m 2 ([/)) 
4deU7 ' 



where the numerical factors are incorporated into the denominators to provide convenient normalization. If U is now 
written in the following form: 

U = e ia k 1 Ak 2 = e iQ fci exp{^(cic^cr 2 + c % a x y a\ + c 3 ^o- 2 )}fc 2 , (29) 
we can compute its local invariants as: 

G\ = cos 2 C\ cos 2 C2 cos 2 C3 — sin 2 Cj sin 2 c 2 sin 2 C3 H — sin 2cj sin 2c 2 sin 2C3 

4 (30) 

o o o o o o 

G2 = 4 cos ci cos C2 cos C3 — 4 sin ci sin C2 sin C3 — cos 2ci cos 2c2 cos 2C3. 

Because the local invariants G\ and G2 characterize the non-local properties of unitary operations, we can use 
these two invariants to check whether a pair of two-qubit gates are locally equivalent. The invariants G± and G 2 
are evaluated by taking the matrix representation of a gate in the Bell basis and then using Eqs. (^6|) and (^). For 
example, CNOT and Controlled-Z (referred as C(Z)) possess identical values of the local invariants, given by G\ = 
and G2 = 1. Therefore, they belong to the same local equivalence class. We refer to this class as [CNOT]. On the 
other hand, the local invariants for VSWAP are G\ — i/4 and G 2 — 0. Hence this gate belongs to a different local 



equivalence class that we refer to as [VSWAP]. Note that from Eq. (]2^), since the local invariants are functions of 
eigenvalues of the matrix to, the local equivalence class can alternatively be defined simply via the set of eigenvalues 
of the matrix to. 



B. Geometric representation of two-qubit gates 

Eq. d3(i| ) reveals the relation between the local invariants G\ and G2 and the coordinates [ci, c 2 , C3] of the 3- Torus 
structure of non-local two-qubit gates. From this relation, given a set of coordinates [ci, C2, C3], we can easily compute 
the local invariants for a local equivalence class. Vice versa, from a given pair of values of the local invariants G\ and 
G2, we can also find the points on the 3-Torus that correspond to a given two-qubit operation. In general, we expect 
to find multiple points on the 3-Torus for a given pair G\ and G2 . We now show how this multiple- valued nature can 
be removed by using the Weyl group to construct a geometric representation that allows the symmetry to be reduced. 

To visualize the geometric structure of the two-qubit gates, we first consider a cube with side length 7r in the vector 
space a. This provides an equivalent representation of the points on the 3-Torus, since T 3 = K 3 /Z 3 . Clearly, every 
point in this cube corresponds a local equivalence class. However, different points in the cube may belong to the same 
local equivalence class. For example, both the points [?,?,?] an d [f > X' T"] corres P on d to the gate VSWAP. 

We use the theory of the Weyl group to reduce this symmetry in the cube. From the Lie group representation 
theory, the orbits of local gates K acting on SU(4)/SU(2) <£> SU(2) are in one-to-one correspondence with the orbits 
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FIG. 2: Illustration of the tetrahedral representation of non-local two-qubit operations. (A) Divide the cube by the planes 
C1—C3 = 0, C1+C3 = 7r, C2 — C3 = 0, and C2 + C3 = n. (B) One of the six equivalent square pyramids produced from (A). Further 
dividing this pyramid by the planes ci — C2 = and c\ + C2 — n gives (C), the tetrahedron OA1A2A3, with Ai = [w, 0, 0], 
A2 = [■§, f , 0], and A3 = [■§,■§,■§]■ OA1A2A3 is a Weyl chamber, denoted a + , with the exception of points on its base where we 
have an equivalence of LA2A1 with LA2O, where L is the point 0, 0]. Every point in a + corresponds to a local equivalence 
class of two-qubit operations. 



{X ea 


Cl 


- c 2 


= 0}, 


{la 


Cl - 


hc 2 


= 


{lea 


Cl 


- c 3 


= 0}, 


{la 


Cl - 




= ^}, 


{leo 


C2 


- c 3 


= 0}, 


{la 


C2 " 


hC3 


= 7T}. 



of the Weyl group WYG,if) on a [41] . From Proposition |[ the Weyl group W(G,K) is generated by the reflections 
s a as given in Eq. (15). Note that in Eq. (|l5|), the reflections s a are either permutations or permutations with sign 
flips of two entries in [c\, C2, C3]. Therefore, if [c±, c%, C3] is an element in a local equivalence class [U], then [a, Cj, c&], 
[w — Ci,ir — Cj, Cfc], [cj, 7t — Cj, 7T — cj,], and [7r — Cj, Cj, 7r — Ck) are also in [[/], where k) is a permutation of (1,2,3). 
With the meaning clear from the context of the discussion, in the remainder of this paper we shall use the triplet 
[ci, C2, C3] to denote either the corresponding local equivalence class of a two-qubit gate, or simply to refer to a specific 
point on the 3-Torus or cube. 

Since each orbit of the Weyl group W(G, K) on a contains precisely one point in a Weyl chamber, the local 
equivalence classes of two-qubit gates are in one-to-one correspondence with the points of a Weyl chamber. Hence, 
each Weyl chamber contains all the local equivalence classes. Recall that the Weyl chambers are obtained by dividing 
the vector space a by the hyperplanes a(X) = 0, where a £ A p as given in Eq. (|l3]). Therefore, we can obtain the 
Weyl chambers by dividing the cube by the planes 

= 0}, {X E a : ci +c 2 = tt}, 

(31) 

Figure ||(A) shows that after dividing the cube by the planes ci — c 3 = 0, c\ + c 3 = tt, c 2 — c 3 = 0, and c 2 + c 3 = n, 
we obtain six square pyramids. One of these pyramids is shown in Figure ||(B). Further dividing this pyramid by 
the planes c\ — c 2 = and c\ + c 2 = tt, we get a tetrahedron OA!A 2 A 3 such as that shown in Figure |(C). Notice 
that for any point [ci,c 2 ,0] on the base of this tetrahedron, its mirror image with respect to the line LA 2 , which is 
[tt — ci, c 2 , 0], corresponds to the same local equivalence class. Therefore, with the caveat that the basal areas LA 2 Ai 
and LA 2 are identified as equivalent, we finally arrive at the identification of the tetrahedron OAiA 2 A 3 as a Weyl 
chamber, and we denote this a + . There are 24 such Weyl chambers in total, and each of them has the volume 7r 3 /24. 
Note that every point in a + corresponds to a different local equivalence class. Consequently, the Weyl chamber a + 
provides a geometric representation of all the possible two-qubit gates. 

For a given two-qubit gate, it is important to find its coordinates [ci, c 2 , c 3 ] on the 3-Torus and hence in the Weyl 
chamber o + . With this representation in the tetrahedron a + we have removed the multiple- valued nature of the 
coordinates on the 3-Torus and cube and therefore can now take the coordinates [ci,c 2 ,c 3 ] as an alternative set of 
local invariants. They provide a useful geometric representation of local invariants that is easy to visualize and is 
entirely equivalent to G\ and G 2 . They can be used directly to implement the local equivalence class of particular 
prescribed two-qubit gates for a given Hamiltonian. More generally, this alternative set of local invariants helps us to 
gain a better understanding of the local invariants and geometric representation of two-qubit gates. 

It is clear that the local gates K correspond to the points O and Ai in Figures || (A)-(C). We now study several 
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nontrivial examples of non-local gates to determine the corresponding coordinates [01,02,03] in a + . All the other 
points of a particular local equivalence class in the cube can be obtained by applying the Weyl group W(G, K) to the 
corresponding point in a + . Note that for a gate [ci, C2, C3] in a + , its inverse is just — c\, C2, C3]. 

(1) CNOT 

Following the procedure to compute the local invariants described above (see Eqs. (|2^) and (^8|)), we obtain G\ = 
and G 2 = 1 for the two-qubit gate CNOT. Solving Eq. @: 

cos 2 ci cos 2 C2 cos 2 C3 — sin 2 ci sin 2 C2 sin 2 C3 = 0, 

sin 2ci sin 2c2 sin 2c3 = 0, (32) 
— cos 2ci cos 2C2 cos 2C3 = 1 , 

we find that [§,0, 0] is the corresponding point for CNOT in the Weyl chamber a + . This is the point L in Figure ||(C). 

(2) SWAP 

For the gate SWAP, we have G\ = — 1 and G2 = —3. Solving Eq. (|30|), we obtain that the corresponding point for 
SWAP is [f , § , § ], i.e., the point A 3 in Figure |(C). 

(3) VSWAP 

The local invariants for the gate VSWAP are Gi = | and G2 — 0. Solving Eq. (|3(]) for this case, we derive that 
[?,?,?] is the corresponding point in a + . This is the midpoint of OA3 in Figure ||(C). 

(4) Controlled- U 

Suppose U is an arbitrary single qubit unitary operation: 

U = exp(7ii(T a; + + 73*0-2). (33) 



For the Controlled-!/ gate, the local invariants are G\ — cos 2 7 and G2 = 2 cos 2 7 + 1, where 7 = By 
solving Eq. |3(]), we find that [7, 0, 0] is the corresponding point in a + . Hence, all the Controlled-!/ gates correspond 
to the line OL in a+, where L is [CNOT]. 



IV. CHARACTERIZATION OF PERFECT ENTANGLERS 



Entanglement is one of the most striking quantum mechanical features that plays a key role in quantum computation 
and quantum information. It is used in many applications such as teleportation and quantum cryptography [ p"5| . In 
many applications, it is often desired to generate maximal entanglement from some unentangled initial states. The 
non-local two-qubit operations that can generate maximal entanglement are called perfect entanglers. In this section, 
we study the perfect entanglers using the geometric approach established in the preceding sections. We will prove a 
theorem that provides a sufficient and necessary condition for a two-qubit gate to be a perfect entangler. It turns 
out that whether a two-qubit gate can generate maximal entanglement is only determined by its location on the 
3- Torus, or more specifically, in the Weyl chamber a + . We show that in the tetrahedral representation of non-local 
gates summarized in Figure ||(C), all the perfect entanglers constitute a polyhedron with seven faces, whose volume 
is exactly half that of the tetrahedron. This implies that the among all the non-local two-qubit operations, precisely 
half of them are capable of generating maximal entanglement from some initially separable states. 

For a two-qubit state ip, define a quadratic function Ent V> = ip T Pip : where P = — ^c^cr 2 0. It can be shown that 
max^, I Ent "01 = 5, and Ent^ = if and only if ip is an unentangled state. The function Ent thus defines a measure 
of entanglement for a pure state. If | Ent-0| = i, we call ip a maximally entangled state. It can be proved that the 
function Ent is invariant under the local operations. 

Definition 3 (Perfect entangler) A two-qubit gate U is called a perfect entangler if it can produce a maximally 
entangled state from an unentangled one. 

Definition 4 (Convex hull) The convex hull C of N points p\, . . . , pn in K" is given by 

, N N \ 

C = <Y^ d jPj I 6 o ■> f° r aU 3 and °i ' = 1 r • ( 34 ) 
^ 3=1 3=1 ' 

Theorem 1 (Condition for perfect entangler) A two-qubit gate U is a perfect entangler if and only if the convex 
hull of the eigenvalues of m(U) contains zero. 
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This result was first mentioned by Makhlin Q but no proof was given. We provide here a proof and then go on to 
develop a geometrical analysis that provides a quantification of the r elative volume of perfect entanglers in SU (4). 

Proof: From the Cartan decomposition of su(4) in Section |ll B| , any two-qubit gate U G U{A) can be written in 
the following form 

U = e la hAk 2 = e M fciexp{^( Cl a>2 +C2CT i CT 2 + C3(T i (J 2 ) | fc2; (35) 

where ki, k% € SU(2) <g> SU(2). For any arbitrary unentangled state ipo, we have 

Ent Uipo = Ent e la k 1 Ak 2 tp = e i2a Ent Atp, (36) 

where -0 = k 2 tp Q is again an unentangled state. From Eq. (p6|), it is clear that | Ent Utpo\ = | Ent Therefore, £/ 
is a perfect entangler if and only if A is a perfect entangler. Furthermore, we have 

Ent Atp = ip T A T PAip 

= (Q^j) t (Q^AQ) t (Q t PQ)(Q^AQ)(Q^j) (37) 

where Q and F are defined as in Eq. ( |l8| ) and (pl|), respectively. The last equality in Eq. (37) holds since Q T PQ = |/. 
Let cj) — Q'ip. Since tp is an unentangled state, we get Ent-0 = 0. Hence, 

Ent = V T PV = <p T Q T PQ<P = \f<P = \(<£ + <t% + & + <Pl) = 0. (38) 
Since tp'tp — lj we have cp^cp = 1, that is, 

|0l| 2 + |02| 2 + |03| 2 + |04| 2 = l. (39) 

Recall the definition of F from Eq. (pl|): 

F = diag{e l = , e l * , e 1 * ,e l 2 }. (40) 



For simplicity, we denote the eigenvalues of F as {Xk}\ =v Then the eigenvalues of m(U) are just {A 2 };j, =1 . We have 

4 

2"' ' • ' ' ' 2 ' * 2 



Ent^ = i(QV,) T F 2 (QtV,) = \<f> T F*<f>=\Ys&\l. (41) 

fc=i 

If A is a perfect entangler, we have 

\ = I Ent A1>\ = ||0?A? + 2 A 2 + cfeXl + $A2| 

<^(I0 2 a 2 | + |^a2| + |^ A 2| + | 2 a2|) (42) 

= ^(I0 2 | + I02l + I03l + I0 2 |)=^ 

The equality in Eq. (^) holds if and only if there exists a real number £ [0, 2tt] such that 

#A? = |* ir V M , <P\\\ = \<p 2 \ 2 e™, 2 A 2 = |0 3 |V 2e , 2 A 2 = |0 4 |V 2e . (43) 
From Eq. d38|), we obtain 

^ , ^ , ^ , , i^i 2 , i^i 2 , i^r 



A 2 A 2 A 2 A 4 



0. (44) 



Since j- = Afc, it follows that 

|0!| 2 A 2 + |0 2 | 2 A 2 + |0 3 | 2 A 2 + |0 4 | 2 A 2 - 0. (45) 
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From the relation in Eq. (j39|), we conclude that if U is a perfect entangler, the convex hull of the eigenvalues of m(U) 
contains zero. 

Conversely, suppose the convex hull of the eigenvalues of m(U) contains zero, that is, there exist {ak}\ = i C [0, 1] 
such that 



a\X\ + a 2 2 X 2 2 + ajXj + a\X\ = 0, 

a l + a 2 + a 3 + a i = 1- 

Let 

. T 

OL\ Ct2 Ck3 tt4 



Ai Aa A3 X4 



(46) 



(47) 



and ip = Qcj). From Eq. (p8|), we have 



_ , 1 , T . 1 /a? os as 

Entip = -4> T (f>= - -| + -| + 3 + 4 =Q ^ 



2 2 y A-^ A2 Ag A^ 

Hence ^ is an unentangled state. From Eq. j4l|), we derive 

Ent Aip = \fF 2 <i> = l -(a\+al + al+al) = 1 -. (49) 

Therefore, U is a perfect entangler. ■ 

We now derive the conditions under which points [01,02,03] in the Weyl chamber a + are perfect entanglers. We 
begin with two corollaries to Theorem |l|. 

Corollary 1 If [01,02,03] is a perfect entangler, then [it — 01,02,03] and [~ — c%, ? — C2, \ — 03] are both perfect 
entanglers. 

Proof: We know that [01,02,03] and [01,-02,-03] correspond to the same two-qubit gate. Since the 3- Torus has 
the minimum positive period 7T, [— tt + Ci, —02, —03] also belongs to the same local equivalence class. From Eq. ( [44] ) 
and (f45|), if [ci, 02, 03] is a perfect entangler, so is [— ci, — c 2 , —03]. Therefore, [tt — 0%, 02, 03] is a perfect entangler. 

From Theorem [l], U is a perfect entangler if and only if the convex hull of the eigenvalues of m(U) contains zero, 
that is, there exist {otk}i = i C [0, 1] such that 

a 2 e t( Cl -c 2 +c 3 ) + a 2 e i(c 1+ c 2 -c 3 ) + a 2 e -i(c 1+ c 2 +c 3 ) + Q 2y (- Cl + C2 +c 3 ) = ^ ^ 

a\+al + al + al = l. (51) 
Substitute the coordinates of the point [^ — Ci, ^ — 02, § — 03] into Eq. (p0[): 

i{a 2 ie ~ i(ci - C2+C3) + a 2 2 e- l{ci+C2 ~ c ^ + o|e i(ci+C2+C3) + a 2 e-^- Cl+C2+c ^ } = 0. (52) 

Together with Eq. ([H]), it is clear that [f — ci, | — 02, § — 03] is a perfect entangler. ■ 

Corollary 2 For a two-qubit gate U, if its corresponding point in the Weyl chamber a + is [ci, £ — ci, 03], [ci,c± — ^, 03], 
or [ci, C2, \ — 02], U is a perfect entangler. 

Proof: For the gate [ci, \ — c±, 03], the eigenvalues of m(U) are 

r „Hc 1 -c 2 +c 3 ) i(c 1 +c 2 -c 3 ) -i(ci+c 2 +c 3 ) gi(-ci +c 2 +c 3 ) j 

' (53) 

= g _i ( c i+ C2 + C3 ) |g i2 ( c i+ c 3) e i2 (ci+c 2 ) ^ e i2(c 2 +c 3 )| 

The convex hull of the eigenvalues of m(U) always contains the origin, and thus [ci, ^ — ci, C3] is a perfect entangler. 
The other cases can be proved similarly. ■ 

Note that for [ci, ^ — 01,03], picking ci = -| and 03 = 0, we obtain the perfect entangler [CNOT]; picking ci — j 
and 03 = j , we get the perfect entangler [ V SWAP] . 

With these corollaries in hand, we can proceed to derive the conditions under which a general point [01,02,03] on 
the 3- Torus is a perfect entangler. 
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e i2(ci+c k ) 



i2(cj+c k ) 



(A) 



(B) 



FIG. 3: Illustration of the proof of Theorem ^| 



Theorem 2 (Perfect entangler on 3-Torus) Consider a two-qubit gate U and its corresponding representation 
[ci, C2, C3] on the 3-Torus. U is a perfect entangler if and only if one of the following two conditions is satisfied: 



7T IT 
-<Ci + C k <Ci + Cj + -< TT, 

3ir ir 

— < Ci + C k < Ci + Cj + — < 27T, 

2 2 



(54) 



where (i,j,k) is a permutation of (1,2,3). 



Proof: 

M g i2(ci+c 2 ) 



Given the eigenvalues of m(U) in Eq. (|53|), it suffices to study whether the convex hull of 
e i2(ci+c 3 ) e i2(c2+c 3 )| con t a i ns the origin or not. Suppose that one of the conditions in Eq. d54J) is satisfied. 
In this case the points {e l2 ( Cl+C2 \ e l2 ' Cl+C3 \ e l2 ( C2+C3 )} have to be on the unit circle as shown in Figure^] (A) or (B). 
It is clear that the convex hull of these three points contains the origin. From Theorem [l], U is therefore a perfect 
entangler. 

Conversely, suppose that U is a perfect entangler. Then the convex hull of the eigenvalues of m(U) contains the 
origin. If all the three points {e* 2 ( cl+C2 ), e l2 ( c i+ c 3) ; e «2(c 2 +c 3 ) j are on ^ e upper or lower semi-circle, the convex hull of 
{1, e * 2 ( c i+ c 2) j e 22(ci+c 3 )^ gi2(c 2 +c 3 )| (jQgg n0 ^ contain the origin. Therefore, we can always pick one point on the upper 
semi-circle and one point on the lower semi-circle such that one of the two conditions in Eq. (|54|) is satisfied. ■ 

The above analysis shows that whether a two-qubit gate is a perfect entangler or not is only determined by its 
geometric representation [01,02,03] on the 3-Torus. Recall that in Section [II B, we show that the local equivalence 
classes of two-qubit gates are in one-to-one correspondence with the points of the Weyl chamber a + , which can be 
represented by a tetrahedron as shown in Figure |2|(C). We are now ready for the final stage of the procedure, namely 
to identify those points in the tetrahedron that correspond to perfect entanglers. 

Consider a two-qubit gate [ci, C2, C3]. As shown in Figure^, in the tetrahedron OA1A2A3, we have ci > C2 > C3 > 0. 
Hence 2(ci + C2) > 2(ci + C3) > 2(c2 + C3) > 0. As in the proof of Theorem @, consider the convex hull of 



= i2(ci+c 2 ) 



j«2(ci+c3) e i2(c2+c 3 ) j._ can identify the following three cases of the gates that do not provide maximal 



entanglement, and are thus not perfect entanglers: 

• If ci + C2 < -|, that is, all the e l2 ( Cj+Ck > are on the upper semi-circle, the gate is not a perfect entangler. In the 
tetrahedron OA 1 A 2 A 3 , c\ + c-i < \ corresponds to the tetrahedron LQPO. 

• If c-2 + C3 > 5, that is, all the e l2( - Cj+Cfe ) are on the lower semi-circle, the gate is not a perfect entangler either. 
This case corresponds to the tetrahedron NPA2A3. 

• From Theorem ||, we obtain that the gates represented by points in the set {X € a + | 2(c± +C3) > 2(c2 +C3) + 7r} 
are not perfect entanglers. This set is the tetrahedron LMNAi. 

The set ofperfect entanglers can thus be obtained by removing these three tetrahedra from OA1A2A3. This is done 
in Figure where it is thereby evident that the polyhedron LMNPQA 2 is the residual set of perfect entanglers. Here 
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FIG. 4: Polyhedron LMNPQA 2 corresponds to perfect entanglers in the Weyl chamber q + (see Figure g (C)), where L, M, N, 
P, and Q are the midpoints of the line segments AiQ, A1A2, A1A3, A3Q, and A2Q, respectively. P corresponds to the gate 
VSWAP, N to its inverse, and L to the CNOT gate. 



the point P corresponds to the gate %/SWAP, N to its inverse, and L to the CNOT gate. Computing the volume of 
the Weyl chamber OA1A2A3 and of these three polyhedra, we have 

^(OAxAaAa) = ^, 

3 3 ( 55 ) 

y-(LQPO) = ^, y(NPA 2 A 3 ) = ^, 7(LMNA 1 ) = ^. 

Therefore, the volume of the polyhedron LMNPQA 2 is 7r 3 /48, which is half of the volume of OA1A2A3. This implies 
that among all the non-local two-qubit gates, half of them are perfect entanglers. Note that the polyhedron LMNPQA 2 
is symmetric with respect to the plane c\ = 7r/2, which provides a geometric explanation of Corollary [ll The points in 
Corollary || correspond to the triangles LMN, LPQ, and NPA 2 , which are three faces of the set of perfect entanglers. 
Also recall that the line OL represents all the Controlled-?/ gates. Hence CNOT, located at L, is the only Controlled-t/ 
gate that is a perfect entangler. Thus we see that the geometric representation provides an intuitive visual picture 
to understand the non-local properties of two-qubit gates, as well as allowing quantification of the weight of perfect 
entanglers. 



V. PHYSICAL GENERATION OF NON-LOCAL GATES 



We now investigate the universal quantum computation and simulation potential of a given physical Hamiltonian. 
We first study the gates that can be generated by a Hamiltonian directly. Generally speaking, these gates form a one 
dimensional subset on the 3- Torus geometric representation of non-local gates. For any arbitrary two-qubit gate, we 
will explicitly construct a quantum circuit that can simulate it exactly with a guaranteed small number of operations. 
Construction of efficient circuits is especially important in the theoretical design and experimental implementations 
of quantum information processing. We assume only that we can turn on local operations individually. Our starting 
point is thus any arbitrary single qubit operation and a two-body interaction Hamiltonian. The single and two-qubit 
operations may be, for example, a sequence of pulses of an optical field that are suitably tuned and focused on each 
individual qubit. The qubits may be represented by either a solid state system such as a quantum dot in a cavity 
fj~5f , or by a gas phase system such as an optical lattice fi^] . 



A. Non-local operations generated by a given Hamiltonian 



In this subsection we investigate the non-local gates that can be generated by a given Hamiltonian H for a time 
duration t, that is, U(t) — expiHt. Recall that I in Eq. ([?]) is the Lie subalgebra corresponding to K, the Lie 
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subgroup of all the local gates. Therefore, i can be viewed as the local part in su(4), and p as the non-local part. If 
the Hamiltonian H contains the non-local part, that is, iH D p ^ 0, then H can generate non-local gates. 

We first consider a Hamiltonian H for which iH is in the Cartan subalgebra a, and then extend to the general 
case. Assume H = \{ci<J x a 2 + c 2 (JyCr^ + c^a^a 2 ). The local equivalence classes of U(t) form a continuous flow on 
the 3-Torus as time evolves. This provides us a geometric picture to study the properties of the gates generated by a 
given Hamiltonian. To illustrate the ideas, we consider the following examples. 

Example 1 (Exchange Hamiltonians) 

(1) Isotropic (Heisenberg) exchange: Hi = |(ct^ct^ + cr^cr 2 + ov^cr 2 ) 
In this case, the two-qubit gate U(t) generated by the Hamiltonian Hi is 

U(t) = expiHt = expi^(alal + a y <T 2 + a\u 2 z ). (56) 

Hence the Hamiltonian Hi generates the flow [|, |, |] on the 3-Torus. The local invariants can thus be computed 
fromEq. My. 



G 2 (t) = — V ; V = 4(cos 6 - - sin 6 -) - cos 3 1 = 3 cost. 



(57) 



We reduce the symmetry of the flow to the Weyl chamber a + , as shown in Figure m We obtain that for t E 



[2fc7r, 2kix + 7r], the trajectory is [|, &, |]; and for t E [2kir + it, 2(k + 1)tt], the trajectory is [|, 7r — |, it — |]. Therefore, 
the flow generated by the isotropic Hamiltonian Hi evolves along OA3A1, which corresponds to all the local equivalence 
classes that can be generated by Hi. Moreover, it can easily be seen that VSWAP and its inverse are the only two 
perfect entanglers that can be achieved by this Hamiltonian. 

(2) Two-dimensional exchange, i.e., XY Hamiltonian: Hi = \{cr x a 2 + &y&y) 

The Hamiltonian H 2 generates the flow [|, |, 0] for t S [2kir, 2kn + tt], and [|,7T — |, 0] for t G [2/c7r + 7T, 2(k + 1)tt]. 
Hence the trajectory evolves along OA2A1. It is evident that H 2 can generate a set of perfect entanglers that 
corresponds to the line segments QA 2 and A2M in a + . Note that A2M represents exactly the same local equivalence 
classes as QA 2 . The local invariants of U(t) are Gi(t) = cos 4 | and G 2 (t) = 1 + 2 cost. 

(3) One-dimensional exchange, i.e., Ising Hamiltonian: H3 = \<Jy(j 2 

The trajectory generated by the Hamiltonian H$ in a + is [|,0,0], which evolves along the line OA^ Hence the 
gates generated by the Hamiltonian H3 are all the Controlled-?/ gates. As noted above, CNOT, located at L, is the 
only perfect entangler that can be generated by this Hamiltonian. The local invariants of U(t) are Gi(t) = cos 2 | and 
G 2 (f) = 2 + cost. 

For any arbitrary H = ^(cia x a 2 + czCTyCTy + c^a],a 2 ), the trajectory on the 3-Torus is [cit, c 2 t, c^t]. If both ci/c 2 
and C1/C3 are rational, the trajectory generated by the Hamiltonian H forms a loop on the 3-Torus. If either ci/c 2 
or C1/C3 is irrational, the trajectory forms a proper dense subset of 3-Torus. 

Next let us consider the case when iH £ p. Recall that we have p = UkeK Adfc(a). Hence for any arbitrary iH E p, 
there exists a local gate k E SU(2) eg) SU{2) such that 

Ad k (iH) = iH a , (58) 

where H a = ^{cia x a 2 + c 2 a y a 2 + czo\a 2 z ). It follows that 

U(t) = exp{iHt) = exp (kUH a kt) = k^ exp ( ^(cio^cr 2 . + c 2 a y a 2 + c 3 ala 2 z )t)k. (59) 



Therefore, the trajectory of U(t) in the Weyl chamber a + is [cit, C2t, c^t}. Eq. (|58j) also implies that H and H a have 
the same set of eigenvalues. We can thus use this property to derive the triplet [01,02,03] explicitly. The following 
example shows how to find the flow in the Weyl chamber a + for a given Hamiltonian H with iH E p. 

Example 2 (Generalized exchange with cross-terms) 

Consider the generalized anisotropic exchange Hamiltonian H = ^(J. xx a x a 2 + J yy <Jy<J 2 + Jxy&l&y + Jyx&y&t) discussed 
in |38] . The eigenvalues of H are 



2 ^ \J i^Jxx ~t~ Jyy) 2 (Jxy J yx) 2 ; \J i^Jxx Jyy) 2 ~i~ {Jxy Jyx^ : 
\J {Jxx Jyy) 2 ~i~ {Jxy Jyx) 2 : \J {Jxx Jyy) 2 ~i~ {Jxy Jyx) 1 



(60) 
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whereas the eigenvalues of H a are 



-{-ci + c 3 + c 2 , -ci - c 3 - c 2 , ci + c 3 - c 2 , ci - c 3 + c 2 }. 



(61) 



Since if and iJ a have the same set of eigenvalues, by comparing Eqs. fl60| ) and (|6lj) and recalling that c\ > c 2 > c 3 > 0, 
we find 



ci 



C2= 2 
c 3 = 0. 



\J {^Jxx H - Jyy)^ 1 {,Jxy Jyx)^ 'sj i^Jxx Jyy)^ [Jxy H~ Jyx)^ 



(62) 



Therefore, the flow generated by this Hamiltonian in the Weyl chamber a + is [cit, c 2 i,0], which evolves in the plane 
OAxA 2 . 

Now we consider the general case when iff € su(4) and H contains both the local and non-local part. To derive the 
trajectory of U{t) = expiHt on the 3-Torus, we first compute the local invariants of U(t) as in Eqs. (Eq) and (Eq): 



Gi(t) 
G 2 (t) 



tr 2 (m(U(t))) 
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tr 2 (m(U(t))) - tr (m 2 (U(t))) 



(63) 



Then from the relation of the local invariants and Cj, we can obtain the flow [ci(t), c 2 (t), c 3 (i)] on the 3- Torus by 
solving Eq. ©: 



G\ = cos 2 Ci cos 2 c 2 cos 2 c 3 — sin 2 c\ sin 2 c 2 sin 2 C3 + — sin 2ci sin 2c 2 sin 2c 3 

r\ ry r\ ry ry ry 

G 2 = 4 cos ci cos c 2 cos c 3 — 4 sin c\ sin c 2 sin c 3 — cos 2ci cos 2c 2 cos 2c 3 . 



(64) 



Example 3 (Josephson junction charge-coupled qubits) 

For Josephson (charged-coupled) qubits |^7j, elementary two-qubit gates are generated by the Hamiltonian Hj = 
— ii?j(cr^ + a 2 ) + (E 2 / E^a^a 2 . If Ej is tuned to aE^, a € K, the local invariants can be obtained 



Gi 
G 2 



1 



(1 + a 2 ) 
1 



J]2 (a 2 (x 2 + y 2 -l) + x 2 ) 2 , 



1 + a 2 



3a z 



1 - 4y 2 a 2 + %a*x*y l + Ax z - Ax'^a 



where 



x = cosa 2 ELt, y — cos y (a 2 + \)aEi,t. 

By solving Eq. (|6~I|), we find that the flow generated by this Hamiltonian on the 3- Torus is: 

c\(t) = ct'EjJ, — uj(t), 
c 2 (t) = a 2 E L t + uj{t), 
c 3 (t) = 0, 



(65) 



(66) 



(67) 



where to(t) = tan -1 y/ '(1 + a 2 y 2 )/ (a 2 — a 2 y 2 ). Since c 3 = 0, the Hamiltonian Hj can reach only those local equiva- 
lence classes on the base OAiA 2 , as shown in Figure [|. Therefore, the Hamiltonian Hj is not able to generate the 
perfect entangler [ V SWAP] . The trajectory generated by Hj in the Weyl chamber a + is shown in Figure |5[ 

For the Hamiltonian Hj to achieve [CNOT], we need to solve Eq. ([35|) for Gi = and G 2 = 1. After some algebraic 
derivations, we find 



2a 2 



(68) 
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3tt 
4 




tt/4 ~T/2 3tt/4 
a = 1.1991, i = 2.7309 



3rr 
4 





3tt 
4 




O 7T^2 3^74 w 6 ~/4 3^74 TT 

a = 0.5, t = 20 a = 0.5, t = 80 



FIG. 5: The flow generated by H.j in the Weyl chamber o 4 



It follows that 



2a 2 cos ylT 



(2fc + 1)tt 
4a 2 £ L ' 
2fc + 1 



(69) 



1, 



where k e Z. When El = I, numerical solution of these expressions shows that the minimum time solution for the 
Hamiltonian Hj to achieve [CNOT] is obtained for a = 1.1991, and the minimum time is 2.7309. 

Note that when the Hamiltonian H contains only the non-local part, that is, iH G p, the flow generated by the 
Hamiltonian on the 3-Torus has constant velocity. However, when the Hamiltonian H contains both the local and 
non-local part, the velocity of the flow on the 3- Torus is usually time dependent as shown in the above example. 



B. Design of universal quantum circuits 



Now let us consider how to generate any arbitrary two-qubit operation from a given two-body Hamiltonian together 
with local gates. The local gates form the Lie subgroup SU(2) £g> SU(2), which also contains all the single qubit 
operations. We will show that by applying the Hamiltonian at most three times, together with four appropriate 
local gates, we can exactly simulate any arbitrary two-qubit gate, i.e., we can implement any SU(4) operation. 
Consequently, if accessibility of any local gate is assumed, this results in satisfying the universality condition needed 
for quantum computation or simulation in a very efficient manner. 

From the discussion in the preceding subsection, we know that a given Hamiltonian is not able to generate any 
arbitrary two-qubit operation simply by turning it on for a certain time period. Generally, the set of the gates that 
can be generated by a Hamiltonian directly is a one-dimensional subset of the 3- Torus. For example, we know that 
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[%/SWAP] can be directly generated from the isotropic exchange interaction between two physical qubits, whereas 
[CNOT] cannot be obtained in this way (unless encoding into multiple qubits is employed j34|). [CNOT] can 
however be achieved by a circuit consisting of two VSWAP and a local gate |3^] . We shall adopt the approach of 
constructing a quantum circuit that contain both non-local gates generated by a given Hamiltonian and local gates, 
and show that this quantum circuit can simulate any arbitrary non-local two-qubit operation exactly with only a 
small number of operations. 

Consider a Hamiltonian H with iH £ p. The gate generated by this Hamiltonian for a time duration t is U(t) = 
cxpiHt. Consider the following prototype quantum circuit 



fc 




U(ti) 




fci 




U(t 2 ) 




fc 2 



k n — 1 




U{t n ) 







where kj are local gates, and n a given integer. Note that the circuit is to be read from left to right. The matrix 
representation of the above quantum circuit is 



fc„t/(t„)fc„_i ■ • • k 2 U(t 2 )k 1 U(t 1 )k . 



(70) 



We will investigate the non-local gates that can be simulated by this quantum circuit. Recall that for this Hamiltonian 
H, there exists a local gate fc G SU(2) ® SU(2) such that Adk(iH) = iH a , where H a = j(c\a].a 2 + c 2 <jy<jy + c^ala 2 ). 
Hence 



U(t) = expiHt = exp(Ad fc t iH a t). 
Let lj = (fcfcj_i • • • fco)^, the quantum circuit (|70|) can now be described as 



(71) 



exp ( Ad/j iHJi) 



exp ( Ad; 2 iH a t 2 



exp ( Ad tn iH a t r , 



knk l n 



and its matrix representation is 

knkUl exp ( Ad/„ iH a t n ) ■ ■ ■ exp ( Ad/ 2 iH a t 2 ) exp ( Ad ;i iH a ti). 



(72) 



We can then pick kj such that lj are in the Weyl group W(G, K). In that case, we have Adi j iH a tj 6 a. Since o is a 
maximal Abelian subalgebra, the quantum circuit in Eq. ((72) is locally equivalent to 



exp ( Ad;„ iH a t n + Ad h iH a t 2 



Ad h iH a h). 



(73) 



Proposition |2| tells us that the Weyl group W(G,K) is generated by the reflections s a given in Eq. (15). Hence for 
a given s a , where a £ A p , there exists a local gate k a such that for any X £ a, Adh a (X) — s a (X). Following the 
procedure in Lemma 2.4, Ch. VII in {y|, we obtain k a as in the following table: 



a 


s a ([ci,c 2 ,c 3 }) 


fc a 


1(03 - c 2 ) 


[ci,c 3 ,c 2 ] 


exp f (|o-i + 




i(c 2 - ci) 


[C2, Ci, C 3 ] 


expf (\o\ + 




«(ci - c 3 ) 


[c3,C 2 ,Cl] 


expf (fcT* + 


2 a y) 


i(c 2 + c 3 ) 


[ci,-c 3 ,-c 2 ] 


ex Pf(|^ - 


2 x) 


i(ci + c 2 ) 


[—C 2 , — Ci, c 3 ] 


expf (fcT* - 


2 u zJ 


«(ci + c 3 ) 


[-C3,C 2 , -Cl] 


expf (|o^ - 


i u y> 



Recall that the flow generated by expiif Q i on the 3- Torus is [cii, c 2 i, c 3 i]. By choosing some appropriate lj from the 
Weyl group W(G, K), we can steer the flow generated by the Hamiltonian. For example, if we want to change the 
flow from [c\t, c 2 t, c 3 i] into [c\t, — c 3 £, — c 2 t], we can simply apply the reflection k i ( C2+c . i y. 



Adfe., , , iH a t 



1 \ / ■ TT \ / 1 ^ 2\ 



2 V 2 

(ci^cr 2 , - c 3 cr^cr2 - c 2 ala 2 z )t. 



(74) 



The following example exemplifies this idea. 
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FIG. 6: The flow generated by the quantum circuit k x exp(iHit2)kl. exp(iHiti) in the Weyl chamber 



Example 4 (Construction of CNOT from isotropic exchange Hamiltonian) 

Consider the isotropic Hamiltonian Hi = j(crla 2 + v^Vy + Our goal is to simulate [CNOT] by a quantum 

circuit containing local gates and two-qubit gates generated by Hi . As shown in Figure |(| the flow generated by 
U{t) = cxpiHit in the Weyl chamber a + is [|, |, |], which evolves along OA 3 for t £ [0,7r]. In the Weyl chamber a + , 



at certain 



the point L ([^, 0, 0]) corresponds to [CNOT]. We then want to switch the flow from 
time instant so that the flow can reach the point L. In order to do that, we can simply apply the reflections Si( C2+C3 ) 



[* * *i to r- 



and s 



in series. The corresponding local gate is thus 



k x = fc 4 ( C3 -c 2 )K. t ( C2+C3 ) = exp 2 (3°* + 2° x ' 6XP 2^2 C7x ~ 2°*' = 6XP ~2° x ' 



(75) 



and we have 



k x exp(iHit)kl = exp Ad kx iH x t = exp -^a\o\ - a\a 2 y - (i\a^)t. 



(76) 



Now consider the following quantum circuit 



exp(iHiti) 




4 




exp(iHit 2 ) 







The flow generated by this quantum circuit evolves along the line OA3 for t € [0, ti], and then switches into a direction 
parallel to the line PL in the plane OA3A1 for t>t\. The matrix representation of this quantum circuit is 



k x exp(iHit 2 )kl. exp(iHiti) 



/ 1 2 



y y 



t 2 exp 



exp 1 



,ti + t 2 i 



hi 1 2 
~lfy°y 



(J z a z )t\ 



ti-t 2 i 



(77) 



Hence the terminal point of the flow is [ t:L + t2 . 



'-]. If we choose t\ 



and t 2 = the terminal point is 



none other than [^,0,0], and thus the quantum circuit simulates [CNOT]. As shown in Figure |6|, the flow generated 
by this quantum circuit is OPL, which goes along the line OP first, and after hitting the point P, it turns to the point 
L along the line PL. Since P is nothing but vSWAP, and t 2 = t\ = ^, we arrive at the known result that [CNOT] 
can be simulated by a circuit consisting of two %/SWAP and a local gate (32| . 

We now derive the following theorem which asserts that when n = 3 the quantum circuit ( f70| ) can simulate any 
arbitrary non-local two-qubit gate. This theorem provides a geometric approach to construct a quantum circuit to 
simulate any arbitrary two-qubit gate from a two body interaction Hamiltonian. 
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Theorem 3 (Universal quantum circuit) Given a Hamiltonian H with iH £ p, any arbitrary two-qubit gate 
U G SU (4) can be simulated by the following quantum circuit 



k 




exp(iiiii) 




h 




exp(i_ffi 2 ) 




k 2 




exp(iHt 3 ) 







where kj are local gates. 



Proof: From Cartan decomposition of su(4) in Section II B , any arbitrary two-qubit gate U € SU (4) can be written 
in the following form: 



U = ki exp{-( 7l cr^CT^ + 72^cr^ + 730 .i a 2)| fcr 



(78) 



where ki and k r are local gates, and 71, 72, 73 € K. We also know that for any given iH £ p, there exists a local gate 
k such that Adfc(iiJ) = iH a , where H a = j(c\o- x a 2 + C20- y a 2 + c 3 o-\o- 2 ), and c\ > c 2 > c 3 > 0. Therefore, the flow 
generated by iH on the 3- Torus is [cit, C2t, c 3 i\. The matrix representation of the above quantum circuit is 



k 3 exp(ii5/i 3 )fc 2 exp(i_ffi 2 )fci exp(iHti)ko. 



Let 



h = (kk kt)\ 
l 2 = {kktk kt)\ 
h = (kk 2 k 1 k kiy, 

the quantum circuit ( |79| ) can be written as 

k 3 kH\ exp ( Ad; 3 iH a t 3 ) exp ( Ad; 2 iH a t 2 ) exp ( Ad (l iH a ti)k r . 
Choose some appropriate local gates fco, ki, and fc 2 such that 



^2 ^i(c3 — c 2 ) ^i(ci+c 3 ) J 

i(c2 — ci) ^i(c 3 — c 2 ) ^i(ci+c 2 ) ? 



i 2 = k 
h = k 

and let k 3 = kil 3 k. It follows that the quantum circuit (|79j) is now 
ki exp (Ad; 3 iH a t 3 ) exp (Ad; 2 i-ff a ta) exp ( Ad (l iH a ti)k r 
= h exp {-{c 3 a l x a 2 x - c^cr^ - c x a z a z )t 3 ) exp {-{-c 3 ala 2 x - c x a y a v + c 2 ala 2 z )t 2 ) 



(79) 



(80) 



(81) 



(82) 



(83) 



= h exp ((cjii - c 3 t 2 + c 3 t3)2 cr x cr x + ( c 2*i - cit 2 - c 2 t 3 )-o- y o- 2 y + (c 3 h + c 2 t 2 - cit 3 )-alo- 2 z ))k r . 
To simulate the two-qubit gate U in Eq. (f7^), we only need to solve the following equation: 






71 




72 




\73 



(84) 



Since 




ci(ci - c 2 c 3 ) + (ci + c 2 )c 3 + (ci + c 3 )c 2 > 0, 



(85) 
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we can always find a solution for Eq. (|84|). Therefore, the quantum circuit ( |79|) can simulate any arbitrary two-qubit 
gate. ■ 

From the above constructive proof, it is clear that together with four appropriate local gates, we can simulate any 
arbitrary two-qubit gate by turning on a two-body interaction Hamiltonian for at most three times. Also note that in 
the proof, the way to choose the local gates ko, ki, and ki is not unique. There are many different ways to choose the 
local gates and time parameters so as to construct the quantum circuit that achieves the same two-qubit operation. 
We therefore can pick the one that is optimal in terms of some cost index such as time. 



VI. CONCLUSION 



In this paper we have derived a geometric approach to study the properties of non-local two-qubit operations, 
starting from the Cartan decomposition of su(4) and making use of the Weyl group. We first showed that the geometric 
structure of non-local gates is a 3-Torus. By further reducing the symmetry, the geometric representation of non-local 
gates was seen to be conveniently visualized as a tetrahedron. Each point inside this tetrahedron corresponds to a 
different equivalent class of non-local gates. We then investigated the properties of those two-qubit operations that can 
generate maximal entanglement. We provided a proof of the condition of Makhlin for perfect entanglers j|] and then 
derived the corresponding geometric description of these gates within the tetrahedral representation. It was found 
that exactly half of the non-local two-qubit operations result in maximal entanglement, corresponding to a seven-faced 
polyhedron with volume equal to one half of the tetrahedron. Lastly, we investigated the non-local operations that 
can be generated by a given Hamiltonian. We proved that given a two-body interaction Hamiltonian, it is always 
possible to explicitly construct a quantum circuit for exact simulation of any arbitrary non-local two-qubit gate by 
turning on the two-body interaction for at most three times, together with four local gates. This guarantees that 
a highly efficient simulation of non-local gates can be made with any Hamiltonian consisting of arbitrary two-qubit 
interactions and allowing control of single qubit operations. 
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